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Abstract 



Using the description in terms of the Hubbard operators hole and spin Green's 
functions of the two-dimensional t-J model are calculated in an approxima- 
tion which retains the rotation symmetry of the spin susceptibility in the 
paramagnetic state and has no predefined magnetic ordering. In this approx- 
imation, Green's functions are represented by continued fractions which are 
interrupted with the help of the decoupling corrected by the constraint of 
zero site magnetization in the paramagnetic state. Results obtained in this 
approach for an undoped 32x32 lattice (the Heisenberg model) and for one 
hole in a 4x4 lattice are in good agreement with Monte Carlo and exact 
diagonalization data, respectively. In the limit of heavy doping the hole spec- 
trum described by the obtained formulas acquires features of the spectrum of 
weakly correlated excitations. 

PACS numbers: 71.10.Fd, 71.27.+a, 74.25.Ha, 74.25.Jb 
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I. INTRODUCTION 



The two-dimensional t-J model is one of the most frequently used models for the de- 
scription of Cu02 planes of perovskite high-T c superconductors (for a review, see Ref. p. 
Together with the numerical methods - the exact diagonalization of small clusters,!'! Monte 
Carlo simulationsa and the density-matrix renormalization-group technique^ - a number of 
analytical methods, such as the mean-field slave-boson! and spin-wave approximations, was 
used for the investigation of the model. The latter method which is based on the spin- 
wave description of the magnetic excitations was shown to be remarkably accurate in the 
case of small hole concentrations and zero temperature .0 This approach was extended to 
the ranges of moderate hole concentrations and finite temperatures,! in particular with the 
use of the spin-wave approximation mo difiedi for short-range order.t9 The positions, sym- 
metry and size of the pseudogaps in the hole and magnon spectra, values of the magnetic 
susceptibility and spin-lattice relaxation rates obtained in this approach are close to those 
observed in photoemission, spin-lattice relaxation and neutron scattering experiments on 
cuprate perovskites.00 

The apparent shortcomings of the spin- wave approximation of the t-J model are the vi- 
olation of the rotation symmetry of the spin susceptibility components in the paramagnetic 
state, the predefined magnetic ordering in the Neel state which serves as the reference state 
of the approximation, and the neglect of the kinematic interaction. In this paper we try 
to overcome these shortcomings by using the description in terms of Hubbard operators. 
Green's functions constructed from these operators are calculated with the use of the con- 
tinued fraction representations following from the Mori projection procedure^ To interrupt 
these otherwise infinite continued fractions we use decouplings of the higher-order Green's 
functions arising in later stages of this calculation procedure. Following the idea of Ref. [13|, 
a correction parameter is introduced in these decouplings to fulfill the constraint of zero 
site magnetization in the paramagnetic state. In this state the obtained components of the 
spin Green's functions are rotationally invariant. The self-energy equations are similar in 
their form to the equations derived in the modified spin-wave approximation.!^ In the case 
of heavy doping the pole in the hole Green's function corresponds to a weakly correlated 
nearest-neighbor band. To check the validity of the obtained equations in the opposite case 
of light doping we have performed calculations for conditions which allow comparison with 
exact diagonalization and Monte Carlo results. We found good agreement of our results with 



the results of Refs. E,14 for spin correlations in an undoped 32x32 lattice and for the hole 



spectral function of a 4x4 lattice with one hole. To gain a notion of the spectral function 
in larger lattices it was calculated in a 20x20 cluster. 



II. DESCRIPTION OF THE MODEL 

The Hamiltonian of the two-dimensional t-J model reads 

H = tn m al a a ma + - Jnm (<4 + ^^m) + %tl, (1) 

nmcr " nm n 

where a ncr = |n<j)(n0| is the hole annihilation operator, n and m label sites of the square 
lattice, a = ±1 is the spin projection, |n<r) and |n0) are site states corresponding to the 
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absence and presence of a hole on the site. If Hamiltonian ([!]) is obtained from the extended 
Hubbard Hamiltonianjl! these states are linear combinations of the products of the respective 
3d x 2_ y 2 copper and 2p a oxygen orbitals.0 We take into account nearest neighbor interactions 
only, t nm = tEa^,m+a and J nm = JJ2 a ^n,m+a where t and J are hopping and exchange 
constants and the four vectors a connect nearest neighbor sites. The spin-| operators can 
be written in the Dirac notations as s„ = | J2<r clncr) (ner| and s„ = \ncr) (n, — <r|. The 



chemical potential /i is included into Hamiltonian (|I|) to control the hole concentration. 
X n = |nO)(nO|. The term — ^ I] na X n X n+a is frequently included into Hamiltonian <^). For 
problems considered below this term leads to an unessential renormalization of the chemical 
potential and therefore it is omitted. The operators a ncr , s£, and X n are the Hubbard 
operators in the space of states of the t-J model. 

The states \na) and |n0) satisfy the following completeness condition: 



J2 |no-)(noi + |nO)(nO| = I. 



(2) 



Using this condition and the above expression for s z n the constraint of zero site magnetization, 
which has to be fulfilled in the paramagnetic state, can be reduced to the form 



s~ l s +l 



0. 



(3) 



where angular brackets denote averaging over the grand canonical ensemble and the hole 
concentration x = (X n ) in the homogeneous state. It should be noticed that in accord 
with the Mermin- Wagner theorem!^ the long-range antiferromagnetic ordering is destroyed 
for any nonzero temperature in the two-dimensional system. Therefore the fulfillment of 
constraint (|) has to be ensured for the considered states. 

The above operators satisfy the following commutation (anticommutation) relations: 



-i „+i 



^"'n^nm) L n' I 
1 



-(T(2 n(T (5 nm , 



{ana, oL,a'} = ( l ~ ^n^n) ^nmC' + S^nmV^i { fl n ff , «m<7'} = 0, 



(4) 



^ n ' X m 



0, [<,X m ]=0, [ 



ncr^nm' 



Notice that the hole creation and annihilation operators do not satisfy the fermion anticom- 
mutation relations. This is the consequence of the exclusion of doubly occupied site states 
due to the strong on-site repulsion [see Eq. (Q)]. 



III. CONTINUED FRACTION REPRESENTATION OF GREEN'S FUNCTIONS 



To investigate the energy spectrum and magnetic properties we shall calculate the hole 
and spin retarded Green's functions 



G{kt) 
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-i'0(*) ({«]«,(*), ak}), D(kt) = us 



i6(t)(st(t),s* 



(5) 
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where = iV -1 / 2 Y, n exp(— 2kn)a ncr , s| = N~ 1 / 2 Y, n exp(— ikn)s^, iV is the number of sites 
and ako-(t) = exp(zift)ak<j exp(— zift). In the considered states G(kt) does not depend on a. 

To calculate the above Green's functions we use their continued fraction representations 
which can be obtained using the Mori projection operator technique.^ Let us consider the 



inner product 



A ■ 



of the operators A and B which is defined in such a manner that the 



following conditions are fulfilled: i) {a A + bB) ■ 



A-C^ 



B-C^ 



, a and b are 



arbitrary numbers; ii) [A, H] ■ B' 



A ■ 



H, fit 



the inner products defined as fit|\ ^ A, B^ 



m 



A ■ B^ = B-A^ 



and 



(ABt) 



dt e"' 



7] -> +0 



. We notice that 



(6) 



satisfy the above properties. Let us divide the result of the commutation of some operator 
Aq with the Hamiltonian into longitudinal and transversal parts with respect to A . The 
transversal part Ai is determined as an operator the inner product of which with A is equal 
to zero. Thus, 



[Ao,H]=E Ao + A 1 , 
where Eq is determined from the condition A\ ■ A Q = 0, 



(7) 



En 



[A ,H]-Al A -Al 



Given Aq and Eq, the operator A± may be found from Eq. (^). The commutator of A\ with 
the Hamiltonian will contain already three terms, 

[A U H] = E 1 A l + A 2 + F A . 

The coefficients E\ and F and the new operator A 2 are determined with the use of the two 



orthogonality conditions 



Ao-A] 



0, i = 0, 1, 



E x = [A U H]-A\ A x -A\ 



[A U H]-Al Aq-AI = A x -A\ A .Al 



where we have used the properties of the inner product. This procedure can be continued. 
In each step of it the coefficients and a new operator are determined by the conditions of 
the orthogonality of this operator to all operators obtained previously. Using the properties 
of the inner product it can be shown that only the n-th, (n+ l)-th andfn — l)-th operators 
appear in the commutator of the n-th operator with the Hamiltonian,H3 



[A n , H] = E n A n + A n+ i + F n -\A n - 



E n 



[A n , H] ■ A n A n ■ A n 



n-l 



(8) 



An ■ Al, 



A n -i ■ A n _ x 



As can be seen, algorithm (§) is the modification of the Lanczos orthogonalization procedure 
which is well known in computational mathematics (see, e.g., Ref. [18] and references therein) 
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and physics!! In Eq. (||), operators A n play the role of mutually orthogonal wave functions 
or vectors in the usual Lanczos procedure. 

Following MoriE! we introduce the projection operator P n which projects an arbitrary 
operator Q on the operator A n , 

PnQ = Q ■ A^ n A n ■ A ni 
and determine the time evolution of the operator A n by the equation 

d n - x 



% M Ant = - p fe)Ki, H], A n>t=0 = A n , 

at k=0 



(9) 



Due to the projection operators in Eq. (|9|) this time dependence of the operator differs from 
the conventional one except the dependence of A . To underline this difference we use the 
subscript notation for the time dependence in Eq. @. Notice also that in accord with this 
equation A nt remains orthogonal to operators A i: i < n for t > 0. Let us divide A nt into 
two parts, 

A n t = R n (t)A n + A nt , R n (t) = A nt ■ A^ n A n ■ A n 

From this definition it follows that A' nt = (1 — P„)A nt . Equations (H) and @) determine the 
time evolution of this operator, 

i±-A' nt = R n (t)A n+1 + n(l - P k )[A' nt ,H]. 

at k=Q 

Solving this equation we find 

A nt = R n (t)A n - i [ dTR n (r)A n+ltt - T . 
Jo 

This result allows us to obtain the following equation for the functions R n (t): 

i-R n (t) = E n R n (t) - %F n f dTR n (r)R n+1 (t - r). 
dt Jo 

After the Laplace transformation R n {oj) = —i Jq 00 dt exp(iut)R n (t) this equation reads 

Bju) = [u-E n - F^+xHr 1 • (10) 

If the inner product is defined as the average of the commutator (anticommutator) 
of operators, the function Rq(uj) = Rq{uj) A ■ A\ coincides with the Fourier transform, 



A 



A 



J^°oo dtexp(iut) ({A A )) , of the commutator (anticommutator) retarded 



Green's functions of the type of Eq. ([|). If the inner product is defined by Eq. 
function Ro(oo) coincides with Kubo's relaxation function 



the 



((Ao\4 



dte^((A \Al)) t , ((A \Al)) t = 6(t) dt'([A (t'),A }). (11) 



From Eq. (|T0|) for all these functions we obtain the following continued fraction repre- 
sentation: 
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U - En 



uj- E l 



(12) 



Thus the recursive procedure (H) in course of which the coefficients E n and F n of the con- 



tinued fraction 
functions. 



12) are determined allows us to calculate Green's or Kubo's relaxation 



IV. THE SPIN GREEN'S FUNCTION 



The direct application of Eqs. (||) and ([12]) to the spin Green's function D(ku), Eq. (^|), 



in the numerator of the contin- 
ty we consider Kubo's relaxation 
defined in Eq. (pTT|) . In this case the inner product ([]) in the numerator 



meets with difficulties because the inner product ^ s k , si k 
ued fraction ( I2|) is equal to zero. To overcome this difficu 
function ((s k si 

of the respective continued fraction is nonzero. After calculating the relaxation function the 
spin Green's function can be obtained from the relation 



((•• 



'-k 



c)) + (*k, *-k) 



(13) 



where we dropped the subscript uj in the relaxation and Green's functions. 

We postpone the calculation of the numerator (s k , slk) of the continued fraction and 

consider its other coefficients. From definition (^|) we find that E (s k , si k ^ 

= and therefore A\ is the Fourier transform of the operator 



'k> 



'k' 



+i„-i 



is{ = - ^mn (5ln - <$lm) S+ 1 S 



E 

mncr 



tmnCT (<5lm - <$ln) aL a m,r = A( + A* 



(14) 



Here the dot over the operator indicates the time derivative. As can be seen, A\ contains 
contributions from spin and hole components A s and A h . Using this result in calculat- 
ing Ri(uj) which is the Laplace transform of the function (v4 lt ,v4{) we neglect the terms 

(A^ t ,A s ^) and (Al t ,A h ^). This approximation is motivated by vanishing values of these 
correlations obtained with the decoupling. Therefore 



(A xu 4)~{A h {t),A h ^) + (Al t ,A 



(15) 



where we have additionally neglected the difference between A\ t and A h (t) (again due to 
zero values of the respective decoupling). In accord with our estimation the influence of 
terms connected with holes in (s k , si k ) and 
small in comparison with the quantity 



At 



A^)) on the spin Green's function is 



*4. k ^4 k Mj even for moderate hole concentrations. 

Therefore in the forthcoming discussion we neglect these terms in (s k , si k ^ and ((*<4 k |*4 k )) 
and consider the time evolution of operators in these quantities as determined solely by the 
Heisenberg part of Hamiltonian (II]). In this approximation the numerator of the continued 



fraction representing ((Ai 



A 



reads 



6 



IS 



k> *-k 



4JC 1 ( 7k 



(16) 



where 7 k = i£ a exp(ika), C p = ^EkTkCk and C k = £ n exp[zk(n _ m )](s+ 1 s m 1 ). For 

= 0. Thus, breaking off the 



Ei we get Ei fis k , — is z _^ = (« 2 s k , — is z _ 
continued fraction on this step we obtain from Eqs. (|12"D , fll3|) , (|i5|) and (|i6|) 



w U 



.4 



4Jd( 7k -r 



2 - 2w nfkw) 



(17) 



where the polarization operator and the excitation frequency are given by 

1 

2 



n(kcj) 



At 



<)) (4, ^k)" 1 , ^k = 4JC 1 ( 7k - 1) (si sQ L . (18) 

To calculate (s k , sf_ k J in the above formulas we notice that in the considered case A 2 = 
s k - (*Sk» -*sik) ( s k> sl k ) X s k and 

16J 2 C 1 2 ( 7 k - I) 2 







H 


Z S k , 2S„ k 



( s k5 s -k) 



0. 



(19) 



We set the above result equal to zero in conformity with the approximation made above in 
the continued fraction where we dropped all terms containing A 2 and operators of higher 



orders. Equation (|19|) can be used for calculating fs k , si k J if the value of 
known. An analogous equation was obtained in Ref. [20 with another method. 



I 6 k , *i_ k 



IS 



We calculate 



i 2 s k , -is z _ k 



in Eq. fll9|) by decoupling terms in the second derivative of 



m' 5 n a l 



T ■ 5 n' 5 m 



'1 ,5 m' 5 n 



(20) 



ttCnm (1 



6„ 



sf where 



In the decoupling we approximate sfs^-s^ by the value 
C nm = (s^Sj^ 1 ). In the last expression we took into account that in accordance with 
Eq. ([|) C nn = | for x = [let us remind that we neglect the influence of holes on the value 

of fs k ,si k )]. Following Ref. [13] the parameter a is introduced to fulfill the constraint of 
zero site magnetization (0) in the paramagnetic state. Before carrying out the decoupling 
it has to be taken into account that terms of Eq. (|20|) in which the site index of the s z 
operator coincides with the site index of s +1 or s -1 operators cancel each other. To verify 
this statement it is necessary to take into consideration that in these terms the operators s„ 
can be substituted by — |, since for the spin-| case = — | + s^s" 1 and s^ 1 ^ 1 = 0. To 
retain this exact cancellation it has to be taken into account before the decoupling. As the 
result we find 

% S\ = Oi ^lm^ln (C mn Si — C*l m S n ^ + J\ n J mn (C*ln s m — C\ m S n 
mn 

+ J2 JL[(l - a)CU(af - s z n ) + aC ln (s z n - s\ 
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and after the Fourier transformation 



2--z 



I S 



where 



4JC 1 ( 7k -l) 



16J 2 a|d|(l- 7 k) (A + l+ 7k ) 



(21) 



C 2 1-a 

A = T^TT + 



3 
4' 



|Ci| 8a|Ci| 
Combining Eqs. (|T|) and fl2T|) we find 

1 = AJa (A + 1 + 7k ) . 



f4, s z _ 



(22) 



In the absence of holes Eqs. ([171) an d (|21|) are close to the equations for the spin Green's 
function and the excitation frequency obtained for the two-dimensional Heisenberg antifer- 



romagnet in Ref. |21| , |22| with the use of the equations of motion for Green's functions and 
Tserkovnikov's formalism,^ respectively. In these works, somewhat more complicated de- 
couplings were used. These decouplings contain several decoupling parameters of the type 
of a which depend on the site indices in the decoupled average. These additional parameters 
allow one to obtain somewhat better agreement with numeric simulations. However, to fix 
the additional parameters exterior data from numerical simulations or the spin-wave theory 
have to be engaged and the theory ceases to be closed. 

As can be shown by the analogous calculation of the transversal spin Green's function 



-i 



, in the paramagnetic state 



-,+1 



(23) 



Thus, the rotation symmetry of the components of the magnetic susceptibility is retained in 
this approach. This fact can be used for the calculation of parameters Cj, C 2 and a in the 
above formulas. From Eq. ( |TTD simplified for the absence of holes, Eq. ( P3"D and the relation 



(4(tK 



due- iuJt e^n B (u)B(kuj) 



we find 



C k = 4/jdl (1 - 7k ) uj^ 1 coth (i/?afc) 



(24) 



(25) 



where B(\tE) = — 7r _1 Im D(kE) is the spin spectral function, ub{E) = [exp(/3E) — 1] _1 and 
(3 = T _1 is the inverse temperature. Substituting this equation in the definitions of C±, C 2 
and in constraint (|3]) we obtain three equations for the three unknown parameters C*i, C 2 
and a. This problem can be reduced to the optimization problem and solved by the steepest 
descent method. 

To check the validity of the approximations made above we used the obtained formulas for 
calculating spin correlations in an undoped antiferromagnet. In Fig. [I] our results obtained 



S 



in a 32x32 lattice for three temperatures are compared with data of Monte Carlo simulations 



performed for the same lattice in Ref. [14]. As can be seen, the agreement is good. However, 
it should be noted that at elevated temperatures in our approximation the spin correlations 
are systematically overestimated in comparison with the Monte Carlo results. 



As follows from Eq. (|21|), for low temperatures and large crystals the spectrum of ele- 
mentary spin excitations is close to the spectrum of spin waves£3 For an infinite crystal and 
T = we found a = 1.70494 and C 2 = -C x = 0.206734. In this case in Eq. (|2T|) the param- 
eter A = and the excitation frequency vanishes in the two points of the Brillouin zone, 
k = (0,0) and (it, tt) (here and below the intersite distance is taken as the unit of length). 
For any nonzero temperature A becomes finite which generates a gap at the (tt, tt) point. 
It can be shownHUll that the gap leads to the exponential decay of spin correlations with 
distance and the respective correlation length is defined by the magnitude of the gap. Thus, 
in agreement with the Mermin- Wagner theorem^ for a nonzero temperature the long-range 
antiferromagnetic order is destroyed in the considered two-dimensional system. 

Now let us calculate the polarization operator n(kcj), Eq. (pT8[). Using the decoupling 
which is equivalent to the Born approximationi and the relations 

roo roo 

(a^(t)aL) = / due- luJt e^n F (uj)A(kuj), (at a a ka (t)) = / due- iu3t n F (u)A(kuj), (26) 

J —oo J —oo 

we find 

7j- roo 

Imn(kcj) = -J2 /k'k / du 1 [n F (J) - n F (J - uo)} A(k' - k, J - cu)A(k'uj ; ), 

U) j^., J —oo 

(27) 

r °° du' ImllYW) 



fC 

ReU{ku) =V 

J — i 



where A(ku) = —n 1 ImG(kc<j) is the hole spectral function, n F {uj) = [exp((3u) + 1] 1 , 
/k'k = 2tiV _1 / 2 (7k' — 7k'-k)(sk, sik) -1 / 2 and V indicates Cauchy's principal value of the 



integral. We notice that Eq. (|27|) is close in its form to the polarization operator obtained for 
the t-J model in the spin- wave approximationJI0 We cannot directly compare the interaction 
constants, because the definition of the spin Green's function in this paper differs from the 
magnon Green's functions in Refs. |5|JTI]|. However, we notice that the spin-wave interaction 
constant and the respective quantity /k'k^k 1 i n Eq. (p7|) are of the same order of magnitude 
and tend to zero linearly with |k| when |k| — > 0. The spin-wave constant behaves analogously 
near the (tt, tt) point, while the quantity /^'k^k 1 does so only in the case of an infinite crystal 
and zero temperature. 



V. THE HOLE GREEN'S FUNCTION 

Now let us consider the hole Green's function. To use the continued fraction repre- 
sentation (|T^) for the anticommutator Green's function G(kt), Eq. (^), the average of the 
anticommutator of operators has to be taken as the definition of the inner product in the 
recursive procedure (||). From the commutation relations (§) we find for the numerator of 
the continued fraction (la^, a kcr } ) = |(1 + x) = <p and for the time derivative 
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ika = *Im (l - S l ^f) S m + S l a m-a ~ ~ J ^ ( aS L S l + S m) «1 -a + A* a lo- (28) 
m m 

With these results we get 

E = ({ia ka , a^}) ({a ka , aL}) = £ k + 

(29) 

£k = (4t0 + etCi^" 1 - 3JF 1 0" 1 ) 7k , // = // + (AtF 1 - 3 JCi)0 _1 , 

where F x = iV _1 Ek7k-^k and F k = E n ex P[^ k ( n ~ m )] (at®™)- 

The estimation of t and J based on the parameters of the extended Hubbard modelil 
gives J/t lying in the range 0.2 — 0.3. For low hole concentrations we can approximate the 
parameter C\ by its value in an undoped lattice. For T = 0.02t in a 4x4 lattice C\ = 0.2119, 
while in a 20x20 lattice C\ = 0.2068. With these parameters the unrenormalized hole 
dispersion can be estimated as ~ — 0.27t7k in the former case and — 0.47t7k in the latter 
case. Thus the first approximation of the recursive procedure describes a band which is much 
narrower than the two-dimensional nearest-neighbor band in the absence of correlations 4t 7 k. 
The reason for this is the antiferromagnetic alignment of spins when the hole movement is 
accompanied by the spin flipping. With increasing the hole concentration G\ — > and the 
unrenormalized dispersion tends to its uncorrelated value. 

The hole Green's function reads 

G(kuj) 



oj — Ek — fx 1 — £(ka>) 

(30) 

S(kc^) = (j)- 1 ((Ax A{)) , Ax = za k(T - (e k + //) a k(J , 



where the difference between A u and Ax(t) was neglected. Due to the mentioned small- 
ness of £ k for low hole concentrations and of J in comparison with t, only the term 
jy-1/2 E lm exp(— ikl) ti m (l — sfsfj + sf a m) _ CT may be retained in A x in the calcula- 
tion of ((Ax The terms in Ax which are linear in spin operators produce the following 
contribution to the self-energy: 



32t 2 ^ ,n 
N<j> \f .1.1- 



ff°° a a n F (-Ux) + n B (u 2 ) 2 , , 

// duJxduj 2 — — ( 7k + 7k-k') A(k-k,Lu 1 )B(k,u 2 ). 

JJ-oo LU — UJ\ — UJ 2 + 17] 



Up to the prefactor this expression coincides with the respective term in the hole self-energy 
calculated in the spin-wave approximation.^ 

The term with three spin operators in Ax produces terms in the self-energy which contain 
two- and three-spin Green's functions. To calculate these functions one would have to solve 
the respective self-energy equations which could be derived in the same way as the equations 
in the previous section. However, such program would essentially complicate the calculation 
procedure. One of the possible ways to overcome this difficulty is to use the decoupling in 
the same manner as we applied it in the previous section, this time in the term with three 
spin operators in A\. However, the comparison with the exact diagonalization data shows 
that this approximation does not give satisfactory results. Another way of simplification is 
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suggested by the above observation that the terms with the one-spin Green's function in 
self-energy (^) are similar to the terms obtained in the spin-wave approximation. This gives 
grounds to suppose that the correcting terms containing two- and three-spin Green's func- 
tions can be approximated by the respective terms of the spin- wave approximation modified 
for short-range antiferromagnetic order. Using the results of Refs. |i0| , |2"5| and Eqs. (|24|), ( |26| ) 
we find 



16irt 2 f°° 

Im£(W) = / duJ ' Kb(-^') + n F {u - uj')] 

N Y k' J—oo 



1 + 7k' 



A + 1 - 7k' 



x 



(7k-k' + 7k) 



1 - 7k' . ( , w v J 1 + 7k' 

+ sgn(w ) (7 k _ k / - 7 k ) J 



A + 1 + 7k' 



A 



7k' 



(31) 



x A(k-k>-w')5(kV), 

„ s r°° do/ImEfW) 
Re£(W) =V 

J — i 



-oo 7T UJ' — UJ 



For low hole concentrations the spin spectral function in the above equation can be 
substituted by its value in the absence of holes, 



B(ku) 




7k 



[5 (uj - u k ) - 8 (uj + lou)) ■ 



2V a VA + l + 7 k 
With this substitution and for low temperatures Eq. ( |3"T| ) acquires the form 

]7Tt 2 [\CA 



(32) 



Im SfkccO 



N(j) 



a 



E 



(7k-k' + 7k) 



' 1 ~ 7k' 
A + 1 + 7 k / 



+ (7 



k-k' 



7k) 



1 + 7k' 



A + 1 - 7k ' 



x [1 + rag (u;k')] A(k — k', uj — co>k 
+ 



(7k-k' + 7k) \l t 1 : 7k ' (7k-k' - 7k) 4 '" ' k 



A + 1 + 7 k 
x nB(uv)^4(k — k', uj + uj^ 



A + 1 - 7k ' 



(33) 



Excluding the numeric prefactor and some other small details this formula is similar to the 
respective formula of the spin-wave approximation .0 

To check the validity of the approximations made we compare the hole spectral function 
calculated using Eqs. ( |30|) and (^) for the case of one hole in a 4x4 lattice with the available 
exact-diagonalization data obtained in this system.! The left panels in Fig. ^| demonstrate 
the results of the exact diagonalization, the right panels present our calculations. Both 
series of calculations were performed for the same set of parameters: J/t = 0.2, t < 0, T = 
and 77 = O.lt (it] was added to the frequency uj in the denominator of Green's function 
to visualize 5-functions; in our calculations parameters C\ and a were estimated for low 
but finite temperature T = 0.02t). As can be seen from the figure, the spectral functions 
obtained in our calculations are in good agreement with the functions found in the exact 
diagonalization. This agreement is somewhat better than that achieved in the spin-wave 
approximation,0 because in contrast to this approximation our approach takes into account 
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the difference between the spectral functions for wave vectors separated by (tt, tt) [cf. the 
spectra for k = (0, 0) and (n, tt)]. As can be seen, in our approximation the binding energy of 
the quasiparticle peak is underestimated in comparison with the exact-diagonalization result. 
This may be connected with the fact that the considered hole concentration x = 1/16 is not 
low enough and can lead to some deviations from the used spin spectral function fl52|) . 

An example of the one-hole zero-temperature spectra in a larger lattice is given in Fig. || 
We notice that the shapes of the spectra cease to change perceptibly with increasing lattice 
size starting from a 16x16 lattice. Excluding the mentioned difference in spectra with wave 
vectors spaced by (n, n) they are close to those obtained in the spin-wave approximation.^ 

With the use of the Hubbard operators equations similar to Eq. ([!!]) were obtained 
also in Ref. |26| for the two-dimensional t-J model and in Ref. ^7| for a somewhat different 
model of the CuC>2 plane. The interaction constant derived in the former work differs 
from the constant in Eq. (|3lD . The constant of Ref. is not applicable for low hole 
concentration: the spectral functions calculated with it differ essentially from those obtained 
by the exact diagonalizationi and in the spin-wave approximation.B! However, this constant 
can be applicable in the region of heavy doping. 



VI. CONCLUDING REMARKS 



Equation (|3l]) was obtained under the supposition of a small hole concentration. We 
have verified that in this limit Eq. (|31| ) in combination with Eq. ( |30| ) describes the hole 
spectral function in good agreement with the exact diagonalization data. On the other 
hand, with increasing hole concentration spin correlations are weakened, the self-energy 
becomes small and elementary excitations described by the two equations tend to the weakly 
correlated nearest-neighbor band with the dispersion e^, Eq. (p9|), where C\ — > 0. Thus, the 
obtained equations give the correct behavior of the hole spectrum in the two limiting cases. 
Besides, it was demonstrated that Eqs. ( |i~7|) and (^T|) with the parameters determined self- 
consistently give a quantitatively correct description of the spin subsystem in the undoped 
case. Equations of the spin- wave approximation, which are similar to Eqs. (|P7D, (f2~l|) and 
(p7|), describe the rapid weakening of spin correlations with hole doping,H3 as it is necessary 
for the above- discussed transformation of the hole spectrum from light to heavy doping. 
This gives ground to suppose that the obtained equations can provide a qualitatively correct 
interpolation between these two limiting cases. 
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FIGURES 



FIG. 1. The spin correlations C(l) = 4|<afsg)|, 1 = (1,0) calculated for T/J = 0.5, 0.75 and 
1 in this work (open circles) and by the Monte Carlo method in Ref. [14| (filled circles). In both 
calculations a 32x32 lattice without holes was used. 

FIG. 2. The hole spectral function AQslj) for the case of one hole in a 4x4 lattice and param- 
eters J = 0.2t, rj = O.lt, and T = 0. Left panels: exact-diagonalization data from Ref. |||, right 
panels: our calculations. The respective wave vectors are indicated in the upper right corners of 
the panels. 

FIG. 3. The hole spectral function A{k.Lo) for the case of one hole in a 20x20 lattice and 
parameters J = 0.2t, ij = O.Olt, and T = 0. Wave vectors indicated near the curves are selected 
along the symmetry lines (0,0) — (0, n) in (a) and (0,0) — (vr,7r) in (b). 
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